Effect of spatial resolution on the estimates of the coherence length of excitons in 

quantum wells 
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We evaluate the effect of diffraction-limited resolution of the optical system on the estimates of the 
coherence length of two-dimensional excitons deduced from the interferometric study of the exciton 
emission. The results are applied for refining our earlier estimates of the coherence length of a cold 
gas of indirect excitons in coupled quantum wells [S. Yang et ai, Phys. Rev. Lett. 97, 187402(2006)]. 
We show that the apparent coherence length is well approximated by the quadratic sum of the actual 
exciton coherence length and the diffraction correction given by the conventional Abbe limit divided 
by 7r. In practice, accounting for diffraction is necessary only when the coherence length is smaller 
than about one wavelength. The earlier conclusions regarding the strong enhancement of the exciton 
coherence length at low temperatures remain intact. 

PACS numbers: 78.67.-n,73.21.-b,71.35.-y 
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I. INTRODUCTION 

Spatial coherence of a bosonic system is encoded in its 
one-body density matrix 



p(r) = (tfVW + r)) 



(1) 



where ('I') is the particle creation (annihilation) op- 
erator and the averaging is over both the quantum state 
and the position r'. In an isotropic system p(r) depends 
on the absolute distance r = |r| only. The density ma- 
trix decreases with r as a result of scrambling the phases 
of the particles' wavefunctions by scattering and thermal 
fluctuations. When this decrease is faster than 1/r 2 , we 
can define a characteristic decay length of p(r) — the 
coherence length — by the relation 



p(r)rdr 




p(r)dr 



(2) 



Coherence length £ a provides a quantitative information 
about fundamental properties of the system of interest. 
For example, a rapid growth of £ x is anticipated as the 
bosons are cooled down below the temperature of their 
quantum degeneracy^ 2 - In addition, £ x sheds light on 
interactions and disorder in the system. 

A gas of indirect excitons in GaAs coupled quantum 
wells is an example of a solid-state system where this 
rich physics can be studied. A number of basic physi- 
cal parameters of such gases (concentration, exciton life- 
time, etc.) can be controlled to bring them to a quasi- 
equilibrium at very low temperatures 4 The coherence of 
excitons is imprinted on the coherence of the light they 
emit4^SiZi2i This allows one to measure £ x by optical 
methods. In this paper we discuss a particular real-space 
technique, which has proved to work in our experimental 
conditions.— It has enabled us to determine the exciton 
coherence length £ x and confirm its rapid increase as tem- 
perature T drops below a few degrees K. 



Traditionally, real-space measurements of the optical 
coherence are done by two-slit (or two-point) interfer- 
ometry. However, this method becomes inaccurate when 
£ x is smaller than the spatial width of the regions from 
which the light is collected. This is the case in our exper- 
iment where £ x does not exceed a few microns. However, 
our technique circumvents this limitation by using a sin- 
gle pinhole. It works well in the regime £ x < D/M\, 
where D is the pinhole diameter and M\ is the image 
magnification factor. In other words, the smallest mea- 
surable £ x is determined not by D/M\ but by the finite 
spatial resolution of the optical system. In this paper we 
show how this resolution can be taken into account. 

The paper is organized as follows. In Sec. [H] we sum- 
marize the main elements of the experimental technique 
and present our principal results. In Sec. IIIII we review 
the theoretical model used in Ref. In Sec. HVI we refine 
it to incorporate the finite-resolution effects. Discussion 
and conclusions are given in Sec. fVl 



II. RESULTS 
A. Measured exciton coherence length 

In order to introduce the key parameters of the prob- 
lem, we need to briefly review the basic steps of the ex- 
perimental implementation of our method for measuring 
the coherence length^ The exciton PL is collected from 
an area of size D/M\ = 2-10 pm in the middle of one 
of the exciton beadsAS ranging ~ 30 pm across. This is 
done by placing a pinhole of diameter D — 10-50 pm at 
the intermediate image plane of magnification M\ = 5, 
see Fig. QJ,. The light is then passed through a Mach- 
Zehnder (MZ) interferometer with a tunable delay length 
51. The output of the interferometer is further magnified 
by the factor Mi w 1.6 (so that the total magnification 
factor is M — M1M2) and then dispersed with a grat- 
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FIG. 1: (a) Experimental setup. The collection angle of the 
lens 2a = 32°. (b) The interference pattern on the CCD 
for D = 25 fim, 81 = 4.2 mm, and T = 1.6 K. (c) Visibility 
function V(T) measured in Ref. H for D — 50 fim and 81 = 
4.2 mm. 



ing spectrometer, resulting in a periodically modulated 
interference pattern. The intensity / = I(x) of this pat- 
tern is recorded by a CCD (Fig. [lb). Here x is the co- 
ordinate along the CCD image. The visibility contrast, 
V = (I max - / ni in)/(W + imin), is calculated (Fig. [It). 
Finally, a theoretical model that relates V to £ x is used 
to analyze the data and determine £ x as a function of 
temperature. 

The results of this analysis are shown in Fig. O As 
one can see, at T < 4K the exciton coherence length 
grows to a few /xm, which exceeds the thermal de Broglie 
wavelength 
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(3) 



by an order of magnitude (AdB ~ 0.1 /xm at T — 2K). 
Here m — 0.2 is the exciton mass in the quantum well in 
units of the bare electron mass. The inequality £ x AdB 
is anticipated for an exciton system near the superfluid 
transition. Possibility of such a transition in systems 
with spatially separated electrons and holes has been put 
forward in Refs.Hlil and [3. The conventional estimate 
of the transition temperature ^ 15 ' 16 Tbkt ~ (^ 2 / m ){n/g) 
gives a few degrees K for the exciton concentrations^ 
n/g ~ 10 10 cm~ 2 (g = 4 is the spin degeneracy). 

In Ref. i we used an approximation of geometrical op- 
tics for describing the light collection in the apparatus. 
This is justified in the most interesting regime of low T 
where £ x is large. On the other hand, at the upper end 
of the temperature range shown in Fig. [2] the estimated 
coherence length £ x was comparable to the diffraction- 
limited resolution of the optical system, e.g., the Abbe 
limiti^ 
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An 



Here NA = sin a is the numerical aperture and Ao is the 
wavelength. In our experiment, NA = sin 16° « 0.3 and 
A w 800 nm. 

In this paper we take diffraction into account and ob- 
tain refined estimates of £ x , which are shown in Fig. [2] 
by the triangles. The difference between previous and 
current estimates is insignificant for all but a few data 
points at the boundary of the experimental resolution. 
Therefore, the case for a rapid and strong onset of the 
spontaneous coherence of the exciton gas below the tem- 
perature of a few degrees K remains intact. 



B. Relation between exciton and optical coherence 
lengths 

It is well understood that the Abbe limit is not a 
"hard" limit but simply a characteristic measure of the 
optical resolution. In fact, another commonly used for- 
mulas due to Rayleigh- '-! differ from Eq. ([4]) by numer- 
ical factors. Roughly speaking, our theoretical model 
enables us to determine which numerical factor is appro- 
priate for our method of measuring £ x . More precisely, 
our main result is as follows. Under certain assumptions, 
^ x is related to the optical coherence length £ by 
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FIG. 2: Main panel: deduced from the fits of V(T) to 
the theoretical curves in Fig. [3] The triangles (squares) are 
the estimates with (without) taking into account the spatial 
resolution of the experimental setup. The dashed lines in- 
dicate boundaries of experimental accuracy. Inset: previous 
estimates of the coherence length.— 
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FIG. 3: Visibility of the interference fringes vs. £ x for pa- 
rameters SI — 4.2 mm, D = 50 pm, Mi = 5, M2 = 1.6, 
Q^ 1 — 0.42 /im. The solid line is the current theory; the 
dashed line is from Ref. 0; the dotted line i s obtained from 
the dashed one by the replacement £ x — > + Q~ 2 - The 
inset shows function V(£ x ) over a larger range of £ x , with the 
dashed line indicating the asymptotic value V(£ x = 00). 



Here, in analogy to Eq. j2]), £ is defined by 

£ = ^ j 9(0, R)RdR I / I / 5(0, -R)<iH , (6) 



where 

g(t, R) = (E(t' + t,R' + R)E(t', R'))/{E 2 {t', R)) (7) 

is the coherence functional of the PL signal E(t, R) emit- 
ted by excitons and collected by the described system. In 
writing this formula we assumed, for convenience, that 
the second magnification (M2) of the image occurs before 
the MZ interferometer, in which case R and R = |R| are, 
respectively, the two-dimensional and the radial coordi- 
nates in the plane of the fully magnified image. 

Equation ([5]) is natural because an experimental mea- 
surement of any length is affected by the spatial resolu- 
tion limit. However, Eq. ([5]) specifically indicates that 
for £ x measured using the setup depicted in Fig. [2^, the 
Abbe limit must be divided by n. This makes its effect 
quantitatively smaller than one would naively think. 

In addition to the limitation from below, £ x > 1/Q, 
the accuracy of the present method is restricted from 
above. When the coherence length exceeds the size of the 
studied region of the sample D/M\, the dependence of V 
on £ x should saturate, see the inset of Fig. [3J This may 
become important at low enough temperatures. The two 
limitations are indicated by the dashed lines in Fig. O 



III. GEOMETRICAL OPTICS APPROACH 

In order to construct the model of the described above 
measurement scheme, we need to know the functional 



form of p(r). Unfortunately, at present there is no 
comprehensive theoretical framework that provides that. 
This is because p(r) is affected by many factors, includ- 
ing thermal broadening and a variety of scattering mech- 
anisms, see more in Sec. [V] On the other hand, the 
present state of experiment 9 - does not allow us to extract 
reliably anything more than a characteristic decay length 
of function p(r). Therefore, we take a phenomenologi- 
cal approach and consider two simple approximations for 
p(r). The first one, used in Ref. [9|, is an exponential 



p(r) = p(0) exp 



(8) 



This ansatz should be reasonable over at least some range 
of r determined by the interplay of the disorder-limited 
mean-free path (~ 1 pm in high-mobility GaAs struc- 
tures), thermal wavelength AdB ~ 0.1pm, and possibly, 
some others. Most importantly, Eq. {8j provides a conve- 
nient starting point because it contains a single charac- 
teristic length £r. Later in Sec. IIVI we will also consider 
a Gaussian form, 



(9) 



p(r) = p(0)exp ( -—^ 



for comparison. 

The crucial point is the relation between p(r) and g(R). 
(Here and below we drop t in the argument of g because 
the time-dependence is not relevant for the calculation.) 
If the experimental apparatus can be described by geo- 
metrical optics, then the only difference between the two 
functions is the linear magnification M and rescaling by 
a constant factor. In this case, Eq. |8j entails 



g(R)=g(0)exp(-R/M^ x ) . 



(10) 



In turn, Eq. © gives £ = £ K , as expected. Until the 
very end of this section we will use these two lengths 
interchangeably 

We have shown previously 9 that the interference visi- 
bility contrast V is related to g as follows: 



dz 



V = 6{1-A)V , 



sin[F(l - A)z] sin[FA(l - z)]g{zD s ) 



V 



f dz 

FA / — sin(Fz)(l - z)g{zD s 



F 



AND S 



D s = M 2 D, A 



SI 



(11) 



where 9(z) is the step-functionjSE A is the linear dis- 
persion of the spectrometer, and N is the number of 
grooves in the diffraction grating. (A = 1.55nm/mm 
and N = 1.5 x 10 4 in Ref. i.) 

For g(r) given by Eq. (|10[) it is straightforward to com- 
pute the integrals in Eq. (| 1 1 [) : however, in general it has 
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to be done numerically. For short coherence lengths, 
£ <C Aq/ ANM, D s /M, one can also derive the analyt- 
ical formula- 



ic 



1 - 



5? 



sin 7rX 



where 



ttX 



Sin = 



\2 



(12) 



(13) 



This equation can be obtained by expanding the sin- 
factors in the integrals to the order 0[z) and extending 
their integration limits to infinity. 

At 51 = SIq and for small enough £, Eq. (112p yields 



1 - 



A 



AND. J D 



M 



(14) 



Thus, as such 51 the visibility contrast V vanishes unless 
£ is nonzero. Working with 51 w 51q ensures the highest 
sensitivity to £. In our experiment^ we extracted £ at 
51 = 4.2 mm, which is close but not exactly equal to 
51q = 5.2 mm. Therefore, we computed V using the full 
formula (fTTj) . The results are shown by the dashed line 
in Fig. [3] We fitted this theoretical curve V(£) to the 
experimentally measured V(T) (Fig. [Ip) using £ = £(T) 
as an adjustable parameter. In this manner we obtain 
the graph shown by the squares in the main panel of 
Fig. [2 We see that the exciton coherence length exhibits 
a dramatic enhancement at T < 4K. On the other hand, 
at T ~ 4K this approach gives £ x — £, ~ Ao- One can 
anticipate that the geometrical optics becomes inaccurate 
at such small £, so that £ and £ K are in fact different. This 
question is studied in the next section. 



IV. DIFFRACTION EFFECTS 

The conventional theory 2 ^, of the image formation in 
optical instruments laid down by Abb©i£ in 1873 predicts 
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FIG. 4: Optical coherence function g — g(R). The solid line 
is computed for the following set of parameters: £ x = 1-5 /im, 
M = 8, Q' 1 = 0.42 fim, and D s = 80 mn. The dashed line is 
for geometrical optics, Q" 1 = 0. 



that a point source imaged by a lens with a magnification 
M creates the diffraction spot 



E{R) cx 



k<Q 



d 2 k 



^^MQ^ (15) 



in the image plane. Here R is the radial distance, Q is 
given by Eq. ©, and J\(z) is the Bessel function. The 
field distribution (fT5)) is known as the Airy diffraction 



pattern.— The physical meaning of Q is the largest tan- 
gential wavenumber admitted by the lens. Accordingly, 
the diffraction can be alternatively viewed as a low-pass 
filtering of the incoming light by the lens^ 

The Airy pattern plays the role of the response func- 
tion of the lens. Its finite spread in R imposes the limit on 
the achievable optical resolution ~ Q^ 1 and is the source 
of the difference between the optical and the actual ex- 
citon coherence lengths, see Eq. (0. Indeed, because of 
the diffraction, the coherence function g(R) is not just 
a rescaled copy of p(R/M) but its convolution with the 
Airy pattern. Using tilde to denote the 2D Fourier trans- 
form, we can express this fact as follows: 



g(k) a. 6(Q - Mk)p(Mk) . 



(16) 



Note that p(k) has the physical meaning of the momen- 
tum distribution function for excitons. Computing p(k) 
from Eq. jSJ), we get 



g(k) = const 



- Mk) 



(1 + M 2 £2fc2)3/2 ■ 



(17) 



The constant prefactor in this formula has no effect on 
V. It is convenient to choose it to be 2ttM 2 ^, so that 



Q/M 



g(R) = M 2 C 



J a (kR)kdk 



{l + M 2 £k 



2)3/2 



(18) 



In this case in the limit Q — > oo, we recover Eq. (|10p with 
g(0) = 1. On the other hand, for finite Q, we have 



5(0) 



1 



1 



(19) 



Additionally, at large R, function g(R) acquires the be- 
havior characteristic of the Airy pattern (| 1 5|) : quasiperi- 
odic oscillations with the envelope decaying as i? -3 / 2 , see 
Fig. |U Finally, computing the optical coherence length 
£ according to Eq. © we get Eq. (0 . 

The refined theoretical dependence of V on £ x can now 
be obtained by substituting Eq. (fl8)l into Eq. ifTT]). As 
before, for small £ x analytical formulas P^| - (fT4^) suffice, 
with £ defined by Eq. (JSJ). When this £ becomes com- 
parable to D s /M, numerical evaluation of Eqs. pip and 
(p~8j) is necessary. The representative results are shown 
by the solid line in Fig. [3) For comparison, two other 
curves are included. The dashed line is the geometrical 
optics approximation, £ — ^ x of Sec. IIII1 The dotted line 
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is the result of correcting the latter according to Eq. ([5]) 
and using Q appropriate for our experiment. As one can 
see, at small £ x , the effect of the diffraction-limited res- 
olution of the optical system is indeed accounted for by 
Eq. ([5]). At large £ x , the correction becomes small and 
all the curves are very close to each other. 

It is instructive to examine how our conclusions so far 
depend on the model assumption (jHJ) about function p(r) . 
To this end we consider next the Gaussian ansatz ([9]), 
which is similar to Maxwcll-Boltzmann distribution ex- 
cept the coefficient in the exponential factor is adjusted 
to satisfy Eq. Let us compute the corresponding £. 
Using Eq. (fTo) , we can rewrite Eq. (O as 



£ = P(fi) I (j Q Q dk~p(k) 
Substituting here 

/ 5(fc)=p(0)exp (~k 2 g) 
which follows from Eq. ©, we get 



crf[(0F/2)Q&] 



(20) 



(21) 



(22) 



where erf (z) is the error function. This formula replaces 
Eq. ([5|). Interestingly, it implies that for the same £ x 
and Q, the effect of the finite resolution in the case of a 
Gaussian decay is always smaller than for the exponential 
one. The direct numerical evaluation of Eq. Ijlip with the 
Gaussian profile © confirms this expectation, see Fig. [5] 
Thus, we again conclude that the diffraction correction 
is important for £ < Ao, but it is very small in the most 
interesting region f > Ao- 

Notice that function V{£ x ) increases somewhat faster 
with £ x for the Gaussian case compared to the exponen- 
tial one, cf. Figs. [3] vs. [5j Therefore, had we adopted the 
Gaussian ansatz ([9]) , the deduced values of £ x (T) would 
have been somewhat smaller than those plotted in Fig. [5] 
This is to be expected: if the exact functional form of p(r) 
is unknown, its characteristic decay length can be deter- 
mined only up to a numerical coefficient of the order of 
unity. 



V. DISCUSSION 

The main purpose of the present work is refinement of 
the optical method for determining the exciton coherence 
length £ x . Unlike previously proposed schemes which 
involve angle-resolved photoluminescence, our technique 
is based on real-space interferometry. 

We showed that in order to obtain an accurate estimate 
of (, x , the optical coherence length £ of the exciton emis- 
sion should be corrected because of the diffraction-limited 
spatial resolution of the experimental apparatus. How- 
ever the correction is insignificant as long as £ x is larger 
than about one wavelength and the numerical aperture 
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FIG. 5: Visibility of the interference fringes vs. £ x for the 
Gaussian coherence function ([9]) and the parameters of Fig. [3] 
The thick solid line includes the diffraction correction; the 
dashed line is for the geometrical optics; the dotted line is 
obtained from the dashed one by the replacement of £ x by 
£ according to Eq. (|22|) . The thin solid line is V(£ x ) for the 
exponential ansatz, replotted from Fig.0to facilitate the com- 
parison. 



NA of the experimental setup is not too small. The cor- 
rection does grow as NA decreases, and so reduced NA 
should be avoided. 

It is well know n 18 i 19 i 21 that limitation of the spatial 
resolution due to diffraction is equivalent to that due to 
restriction on tangential wavenumbers k admitted by the 
lens collecting the signal. This fc-filtering effect has been 
considered in Ref. [2J in application to the measurements 
of the exciton coherence length. For the collection an- 
gle a — 16° in our experiments, 9 - the results presented 
in Fig. 3c of Ref. [13 give the correlation length due to 
the /c-filtering effect £ 7 « 1 p.m. (This length plays the 
role similar to that of Q' 1 — 0.42 /im in our formal- 
ism.) For the considered p(r) this correction enters either 
through the quadratic sum, Eq. ([5]). or the error function, 
Eq. (|2"2"|) . As a result, the estimation of the fc-filtering ef- 
fect per Ref. [U gives only a small (~ 10%) correction, 

e -g-j £ — \J^ 2 ~ £7 to the optical coherence length £ mea- 
sured in Ref. d at low T. Therefore, it cannot explain 
the observed large enhancement of the coherence length 
at T < 4 K. Our calculations indicate that the correction 
is even smaller. 

The discussion of physics that is responsible for the 
observed rapid change in £ x at low temperatures is how- 
ever beyond the score of this paper. As a final word, we 
would like to offer only the following minimal remarks on 
this matter. 

The density matrix p(r), from which £ x is defined 
through Eq. @, is influenced by a number of factors, 
including Bose statistics, interactions, and scattering. 
The effect of the first two has been studied extensively, 
albeit for simplified models of interaction. According 
to present understanding;^ the long-distance behavior 
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of function p(r) is qualitatively different above and be- 
low the Berezinskii-Kosterlitz-Thouless (BKT) transition 
temperature Tbkt- At T 3> Tbkt, where p(k) obeys the 
classical Boltzmann statistics, p(r) decays as a Gaussian, 
Eq. ([9]), with the coherence length 

£x=A dB Ar. (23) 

[Our estimates of £ x at T < 4K exceed Xab/k by an 
order of magnitude, suggesting that Eq. (f2"o) does not 
apply at such temperatures.] At T < Tbkt, the eventual 
asymptotic decay of the density matrix becomes alge- 
braic, p(r) oc r~~ u with a temperature-dependent expo- 
nent v{T). The behavior of p(r) at intermediate temper- 
atures and/or distances is more complicated. In general, 
it can be computed only numerically, e.g., by quantum 
Monte-Carlo method^ 

Some of the other mechanisms of limiting the co- 
herence length £ x , such as exciton recombination and 
exciton-phonon scattering are too weak to significantly 
affect the large magnitude of observed £ x at low tem- 
peratures^ However, scattering by impurities and de- 



fects should be seriously considered. It can substan- 
tially modify the functional form of p(r) compared to 
the disorder-free case. Indeed, weak disorder typically 
leads to an exponential decay of the correlation func- 
tions on the scale of the mean-free path, which in fact 
inspired our ansatz ([5]). As temperature goes down, the 
strength of the disorder decreases because excitons can 
screen it more efficiently j 25 i 26 i 30 This should increase both 
the mean-free path and the exciton coherence length. 

The comprehensive theoretical calculation of the ex- 
citon coherence length that would take into account all 
relevant thermal, interaction, and disorder screening ef- 
fects is yet unavailable. 
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